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Abstract. - We study the electron thermal transport in granular metals at large tunnel 
conductance between the grains, (;t 3> 1 and not too low a temperature T > grS, where 5 
is the mean energy level spacing for a single grain. Taking into account the electron-electron 
interaction effects we calculate the thermal conductivity and show that the Wiedemann- Franz 
law is violated for granular metals. We find that interaction effects suppress the thermal 
conductivity less than the electrical conductivity. 



The electrical conductivity a and the thermal conductivity k of an ordinary metal are 
related to each other via a universal relation [1] known as the Wiedemann-Franz law: 



where Lq = 7r^/3e^ is the Lorentz number, T is the temperature and e is the electron charge 
(we use units such that ks = h = 1). The Wiedemann-Franz (WF) law holds as long as 
scattering processes are (quasi-) elastic i.e. only weakly energy-dependent. 

Equation J^l is a direct consequence of Fermi liquid theory. This theory describes the 
properties of a weakly disordered, interacting electron gas, provided that screening renders 
the Coulomb interactions sufRciently weak and short-ranged. Under these conditions, the 
low-lying excitations of the metal are non-interacting fermionic quasi-particles that carry 
both charge and energy, and as a consequence the WF law is valid [2]. In other words, no 
additional information is obtained from a measurement of the thermal conductivity that is 
not already present in the electrical conductivity as long as interaction effects result only in 
the renormalization of electron spectral parameters. 

By now it is well known that the simultaneous action of disorder and interactions changes 
the behavior of metals at low energies [3] . Interference of diffusively scattered electron waves 
weakens the screening properties of the electron gas, thereby leading to an energy-dependent 
increase of the interaction strength. Therefore, in contrast to an ordinary Fermi-liquid, both 
transport and thermodynamic quantities acquire an additional non-trivial energy-dependence. 
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An example is the so-called zero-bias anomaly: close to the Fermi level Coulomb effects 
suppress the tunneling density of states. Another example is the temperature-dependent 
interaction correction to the Drude conductivity. 

Since the low-energy properties of disordered metals with interactions deviate from those 
of an ordinary Fermi-liquid, one is, in view of the above, naturally lead to the conclusion that 
the Wiedemann-Franz law (or deviations from it) may serve as an important probe uncovering 
the underlying microscopic nature of disordered metallic systems. Although this problem has 
been addressed by several authors [4-8] over the past years, it has not been solved completely 
as yet. Ref. [4] finds that the WF law holds, whereas in Refs. [5-8] deviations from the WF 
relation are found (although these papers do not agree on the precise form of the deviations). 
We refer the reader to Ref. [7, 8] for a detailed account of the present status of the field. 

In this Letter we investigate this issue from a somewhat different perspective consider- 
ing thermal transport in granular metals. A granular metal consists of small metallic grains 
coupled via tunnel junctions [9]. The electron tunnelling between grains induces random- 
ness while the small electrical capacitance associated with the tunnel-coupled grains induces 
Coulomb charging effects. The interest in the physics of granular metals is two-fold. First, 
their properties are generic for a wealth of disordered strongly correlated systems. On the 
other hand, as the interaction strength and degree of disorder can be controlled by properly 
choosing the parameters of the tunnel junctions and the grains, these systems offer a unique 
experimentally tunable disordered interacting system. 

Depending on the value of the dimensionless tunnel conductance gx between the grains, 
we can distinguish a low conductivity regime, corresponding to gx ^ ^ and a metallic, highly 
conducting regime, where (/t 3> 1. Electric transport in the limit ^ 1 has been studied 
extensively [10, 11] since the early work by Mott [12]. More recent work deals with the 
electronic properties of granular systems in the metallic regime [13-16]. 

In what follows we focus on the electron thermal transport in the metallic regime, where 
the dimensionless tunnel conductance is large, gr ^ 1- We show that in the absence of 
interactions the Wiedemann-Franz law holds and the thermal conductivity is given by 

kq = LoaoT, (2) 

where cto — 2e^gT0,^~'^ is the classical Drude conductivity for a granular metal (including 
spin), a is the size of a single grain, and d is the dimensionality of the system. 

We find further that interactions violate the Wiedemann-Franz law in granular metals. 
Specifically, our result for the thermal conductivity can be conveniently formulated in terms 
of the quantity 

Sk = K — LqT a , (3) 

that reflects deviations from the Wiedemann-Franz law. Including interaction effects, Sk in 
Eq. |j3Jl is given by 

{jT/a d = 3, 

where Ec is the charging energy of an isolated grain and 7 « 0.355 is a numerical constant. 
Our main result @ is valid at not too low temperatures, T > gxS, where 5 is the mean 
level spacing of the grains. At these relatively high temperatures [15], the electronic motion 
is coherent within the grains, but coherence does not extend to scales larger than the size of 
a single grain. Under these conditions, the electric conductivity a in the right hand side of 
Eq. ||3Jl is given by the expression obtained in Refs. [15, 16] 



a = ao[l - {l/2TTdgT)H9TEc/T)], 



(5) 
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where the electric conductivity ctq was defined below Eq. Contrary to the electric con- 
ductivity of granular metals, a given by Eq. (jSJ, the thermal conductivity 5k, in Eq. Q is 
sensitive to the long-ranged character of Coulomb interaction and therefore depends on the 
dimensionality of the sample, even at temperatures T > grS. 

One can see from Eqs. (PJ and Q that the thermal conductivity increases with temperature 
in both two and three dimensions. From Eq. Q it follows that for granular films both Sk and 
the correction to the electric conductivity (second term in the right hand side of Eq. (jSJ) have 
the same logarithmic behavior such that in 2d the thermal conductivity can be written as 



(^-2) ^, 

K2d = kq Tin 

6 



(^) , (6) 



where kq is defined by Eq. ©. 

Having summarized the main results, we now turn to the quantitative description of our 
model and derivation of Eq. We consider a d-dimensional array of metallic grains. The 
motion of electrons inside the grains is diffusive and they can tunnel between grains. Inter- 
action effects stem from the long-range part of the Coulomb interaction, which is just the 
classical electrostatic charging energy of the grains. Such an array of weakly coupled metallic 
grains is described by the Hamiltonian H = ^ ^ Hi , where the sum is taken over all grains in 
the system and 

g2 

k j 



"'aLa.... (7) 



The first term in the right hand side of Eq. Q describes free electrons in the i-th isolated 
disordered grain, a|j,(ai^fc) are the creation (annihilation) operators for an electron in the state 
k and = /2'm — /i with ^ being the chemical potential. The second term describes the 
charging energy with the capacitance matrix Cij and the electron number operator fii in the 
i-th grain. The last term represents the tunneling between the grains, where t^J is the tunnel 
matrix element corresponding to the points of contact of i-th and j-th grains. 

The thermal conductivity of a granular metal is given by the analytical continuation of 
the Matsubara thermal current-current correlator. For granular metals the thermal current 
operator J^j between the i-th and the j-th grain can be obtained as follows. The energy 
content of a given grain i changes as a function of time, such that 

dH,/dt^i[H,H,]. (8) 

Energy conservation requires that this energy flows to the other grains in the system. From a 
straightforward calculation of the commutator in Eq. we obtain J^^ , the amount of energy 
flowing per unit time from grain i to a neighboring grain j. It can be written as 



where J^}^ is the heat current in the absence of electron-electron interaction effects 
''J 



J^^ = \Yj^v+Q^<r^3,, + h.c^ (9b) 
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Fig. 1 - Diagrams representing: (a) the diffusion propagator of granular metals 7?(n„,q) in Eq. 
The tunneling vertices are described by the circles; (b) interaction vertex renormalized by impuri- 
ties and intergranular scattering; (c) the effective screened Coulomb interaction for granular metals 
l/(n„,q) in Eq.im 



while the second term, j^J"^ on the right hand side of Eq. (|9a|l appears due to Coulomb 
interaction effects 



j{2) _ _1 



(9c) 



In equation H9c(l the operator is the electric current which is given by the foUowing expres- 



sion 



(9d) 



For large tunnel conductance, the Matsubara thermal current-current correlator can be 
analyzed by virtue of the perturbation expansion in I/qt, using the diagrammatic technique 
developed in Refs. [13, 14] which we briefly outline here. The average electron Green function 
has the same form as that of a homogenously disordered metal with self energy Y,{uin) = 
— (i/2T) sgncjn, where ojn is a fcrmionic Matsubara frequency. The elastic mean free time r is 
mainly determined by the processes of scattering by impurities inside the grain. Intergranular 
scattering processes (tunneling) are included assuming that tunneling elements are random 
Gaussian variables correlated according to 



VlP2 ^f3P4^ _ 



(10) 



Here i and j are the indices of the neighboring grains and t is the tunnel amplitude, which 
is related to the dimensionless tunnel conductance, gr = 27rt^/i5^. Inclusion of intergranular 
scattering processes results in a small renormalization of the mean free time r. The diffusive 
motion inside a single grain is given by the ladder diagram that results in the diffusion prop- 
agator Dq^ ~ T\Qn\ where f2„ is a bosonic Matsubara frequency. The coordinate dependence 
in Dq is neglected since in the regime under consideration all characteristic energies are less 
than the Thouless energy. The complete diffusion propagator is given by the ladder diagrams 
shown in Fig. la resulting in the following expression 



i?"i(r!„,q) -r(|r!„|+£q5). 



(11) 



where eq — "^gr^^i^ — cosqa) with a being the lattice vectors is the function of quasi- 
momentum q that appears due to the tunneling between the grains. The same ladder diagrams 
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(a) (b) (c) 




(d) (e) 



Fig. 2 - Diagrams describing the thermal conductivity of granular metals at temperatures T > grS: 
the diagram (a) corresponds to kq in Eq. Q. Diagrams (b)-(e) represent first order corrections in 
l/gfT to the thermal conductivity of granular metals due to electron-electron interaction. The solid 
lines denote the propagator of electrons, the dashed lines describe effective screened electron-electron 
propagator and the black triangles describe the elastic interaction of electrons with impurities. The 
tunneling vertices are described by the circles. 



describe the renormalized interaction vertex in Fig. lb. This vertex is used to obtain the 
polarization operator, that defines the effective dynamically screened Coulomb interaction for 
granular metals. Fig. Ic, 



\VL.n\+5e 



(12) 



where C(q) is the Fourier transform of the capacitance matrix. 

In the absence of electron-electron interactions, the thermal conductivity is given by the 
diagram (a) in Fig. |21 A direct calculation of this diagram confirms that the thermal conduc- 
tivity is given by Kq in Eq. ^ and thus that the WF law holds. At temperatures T > gxS 
the first order interaction corrections to the thermal conductivity are given by the diagrams 
(b-e) in Fig. 2. We neglect the diagrams that can be obtained from (b-e) by insertion of 
an extra diffusion propagator since their contribution is small in the parameter qtS/T. The 
diagrams (b) and (c) in Fig. 2 stem from contributions of the non-interacting part j'l^^ to the 
thermal current-current correlator. These diagrams are analogous to the ones considered in 
Ref. [16] for the interaction correction da to the conductivity of granular metals. In fact, we 
find that the contribution of diagrams (b) and (c) to 6n in Eq. Q is cancelled exactly by the 
corresponding term LgTSa. Thus, although these diagrams lead to corrections to both a and 
K, they do not lead to a violation of the WF law. Diagrams (d) and (e) contain the interacting 
part jy of the thermal current. There are no analogous diagrams for electric conductivity, 
hence any non-vanishing contribution from (d) and (e) inevitably leads to a violation of the 
WF law. The diagram (d) in Fig. 2 turns out to be zero. Thus, among the diagrams (b-e) in 
Fig. 2 only the diagram (e) contributes to 6k in Eq. Q giving rise to the correction 

27rT2 sinh2(r!/2T) ' ^ ' 

where the sum is taken over the quasi-momentum q and the potential Vq{^l) is given by the 
following expression 

T> (c,\ - [l + cos{q^a^)]Ec{q) , , . 
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Here i?c(q) = e^/2C(q) is the charging energy. For the 3d case we simphfy the right hand 
side of Eq. Ijl3b|) taking the Umit Ec{<i) — > oo- Performing integration over and summation 
over q in Eq. (|13a|l we obtain our final result for the thermal conductivity Q in the 3d case, 
where the numerical coefScient 

_ TTO^ r d^q 1 + cos(g^a^) 
^"67 (27r)3 Ea[l - cos(qa)] ' ' 

For the 2d case one can see that the main contribution to the integral over q comes from 
the region of small quasi-momenta q ^ a^^ where one can use the asymptotic expression 
for the capacitance matrix C(q) — a^g/27r. Performing the integrations in Eq. H13a|) with 
logarithmic accuracy we obtain the final result for the correction to the thermal conductivity 
of 2d granular metals 

So far, we ignored the fact that electron-electron interactions also renormalize the chemical 
potential fi. Generally, this renormalization may affect the thermal conductance: the thermal 
current vertex, Eq. (|9 bp . as well as the electron Green functions depend on /i. To first order 
in the interactions, the renormalization of only leads to corrections to diagram (a) in Fig.|21 
As it can be easily shown, for this diagram the renormalization of the two heat current vertices 
is exactly cancelled by the renormalization of the two electron propagators. Therefore, the 
renormalization of the chemical potential by the interactions does not affect our results in 
lowest order. 

In this paper we only consider the electron contribution to thermal conductivity. As it is 
well known, phonons will provide an independent, additional contribution to thermal trans- 
port. Although it is generally argued that phonons can be ignored at low temperatures, this 
point needs to be considered with care as we will now argue. Consider acoustic phonons that 
propagate with velocity Cph in a three-dimensional impure host, characterized by a phonon 
scattering length Iph- In this simple model, the phonon contribution to the thermal conductiv- 
ity is given by Kph = T^lpu/c^h- We see that Kph kq at a temperature T* ^ gTC^f^/lpha. 
The electron contribution to thermal transport dominates only for T < T*. In particular, if 
the phonon contribution is to be ignored in the temperature range relevant to us, we need to 
impose qtS < T*. It is clearly a matter of system parameters whether or not this condition 
will be satisfied in a typical experiment. Avoiding the phonon contribution in a measurement 
of thermal conductivity thus implies constraints on the choice of materials as well as on sample 
fabrication. 

In conclusion, we have investigated the electron thermal transport in granular metals in 
the limit of large tunneling conductance between the grains and temperatures T > grS. We 
have calculated the thermal conductivity of granular metals and shown that the Wiedemann- 
Franz law in granular metals is violated if electron-electron interaction effects are taken into 
account. Our result can be expressed conveniently through the Lorentz number of the granular 
metal, L — n/aT. In the absence of interactions, the Lorentz number is a universal constant, 
L = Lq ^ 7r^/3e^. In the presence of interactions, we find that L = Lq + SL, where SL = const 
for a three-dimensional granular system, whereas 6L ~ \ii(gTEc /T) for a granular film. In 
both cases, 6L > 0, indicating that interactions suppress electric conductivity more than the 
thermal conductivity. 

Since interacting granular metals behave differently from homogeneous disordered metals, 
we can not directly compare our results with the ones obtained in Refs. [4-8]. Nevertheless, 
the structure of our result for thermal conductivity, as given by Eqs. JSJ and Q, is consistent 
only with that of the results of Ref. [7,8]. Specifically, we find two contributions of opposite 
sign related to interactions. One contribution is the interaction correction to electric conduc- 



I. S. BelOBORODOV et al: THERMAL TRANSPORT IN GRANULAR METALS 



7 



tivity, multiplied by TLq. The other contribution leads to a violation of the WF law and is 
logarithmic in temperature for the 2d case, as in Ref. [7]. Of course the prefactors are different, 
and a complete agreement can be expected only in the low temperature limit, T < grS, which 
is beyond the scope of this Letter. 

Experiments on thermal transport are important since they allow for a direct observation 
of deviations from ordinary Fermi liquid behavior. However, measurements of the electron 
thermal conductivity of a mesoscopic conductor at low energy is still an experimental challenge. 
So far, only two-dimensional electron gas systems were studied [17]; more recent experiments 
aimed at achieving the accuracy that is necessary to observe deviations from the WF law [18]. 
We hope that our results that reveal deviations from Fermi liquid behavior of granular metals 
will stimulate further experiments on thermal transport. 
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